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Abstract 



We consider sequences of random variables whose probability generating functions 
are polynomials all of whose roots lie on the unit circle. The distribution of such random 
variables has only been sporadically studied in the literature. We show that the random 
variables are asymptotically normally distributed if and only if the fourth normalized (by 
the standard deviation) central moment tends to 3, in contrast to the common scenario for 
polynomials with only real roots for which a central limit theorem holds if and only if the 
variance goes unbounded. We also derive a representation theorem for all possible limit 
laws and apply our results to many concrete examples in the literature, ranging from com- 
binatorial structures to numerical analysis, and from probability to analysis of algorithms. 



The close connection between the location of the zeros of a function (or a polynomial) and the 
distribution of its coefficients has long been the subject of extensive study; typical examples 
include the order of an entire function and its zeros in Analysis, and the limit distribution 
of the coefficients of polynomials when all roots are real in Combinatorics, Probability and 
Statistical Physics. We address in this paper the situation when the roots of the sequence of 
probability generating functions all lie on the unit circle. While one may convert the situation 
with only unimodular zeros to that with only real zeros by a suitable change of variables, 
such root-unitary polynomials turn out to have many fascinating properties due mainly to the 
boundedness of all zeros. In particular, we show that the fourth normalized central moments 
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thanks ISM for its hospitality and support. 

^Part of the work of this author was done during his visit at Institute of Statistical Science, Academia Sinica, 
Taipei. He gratefully acknowledges the support provided by the institute. 
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are (asymptotically) always bounded between 1 and 3, the limit distribution being Bernoulli if 
they tend to 1 and Gaussian if they tend to 3. 

Although this class of polynomials does not have a standard name, we will refer to them as, 
following Kedlaya (2008) and for convention, root-unitary polynomials. Other related terms 
include self-inversive (zeros symmetric in the unit circle), reciprocal or self-reciprocal (P(z) = 
z n P(z~ 1 )), uni-modular (all coefficients of modulus one), palindromic (aj = a n _j), etc., when 
P{ z ) = So<j<n a j z ^ * s a polynomial of degree n. 

Unit roots of polynomials play a very special and important role in many scientific and en- 
gineering disciplines, notably in statistics and signal processing where the unit root test decides 
if a time series variable is non-stationary. On the other hand, many nonparametric statistics are 
closely connected to partitions of integers, which lead to generating functions whose roots all 
lie on the unit circle. We will discuss many examples in Sections 4 and 5. Another famous 
example is the Lee- Yang partition function for Ising model, which has stimulated a widespread 
interest in the statistical-physical literature since the 1950's. 

While there is a large literature on polynomials with only real zeros, the distribution of the 
coefficients of root-unitary polynomials has only been sporadically studied; more references 
will be given below. It is well known that for polynomials with nonnegative coefficients whose 
roots are all real, one can decompose the polynomials into products of linear factors, implying 
that the associated random variables are expressible as sums of independent Bernoulli ran- 
dom variables. Thus one obtains a Gaussian limit law for the coefficients if and only if the 
variance tends to infinity; see the survey paper Pitman (1997) for more information and for 
finer estimates. A representative example is the Stirling numbers of the second kind for which 
Harper (1967) showed that the generating polynomials have only real roots 1 ; he also established 
the asymptotic normality of these numbers by proving that the variance tends to infinity. For 
more examples and information on polynomials with only real roots, see Brenti (1994), Pitman 
(1997) and the references therein. See also Haigh (1971), Hayman (1956), Renyi (1972) for 
different extensions. 

Our first main result states that if we restrict the range where the roots of the polynomials 
P n (z) can occur to the unit circle \z\ = 1, then the asymptotic normality of X n defined by the 
coefficients is determined by the limiting behavior of its fourth normalized central moment. 
Throughout this paper, write X* := (X n — E(Jf n ))/ a/V(X„). 

Theorem 1.1. Let {X n } be a sequence of discrete random variables whose probability gen- 
erating functions K(z Xn ) are polynomials of degree n with all roots pj lying on the unit circle 

Pi. ~ 1- 

- (Bounds for the fourth normalized central moment) 

1^E(X:) 4 <3 (n^l). (1) 

- (Asymptotic normality) The sequence of random variables {X*} converges in distribution 

(and with all moments) to the standard normal law ,sY((), 1) if and only if 

E(X*) 4 ^3 (n->oo). (2) 

'The fact that the Stirling polynomials (of the second kind) have only real roots had been known long before 
Harper; see for example d'Ocagne (1887); in addition, Bell (1938) wrote (without providing reference) that all 
results of d'Ocagne's paper were already obtained thirty years before him by a number of English authors. 
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- (Asymptotic Bernoulli distribution) The sequence {X*} converges to Bernoulli random vari- 
able assuming the two values — 1 and 1 with equal probabilities if and only if 

E(X*,) 4 ^1 (n-»oo). (3) 

This theorem implies that Gaussian and Bernoulli distributions are in a certain sense ex- 
tremal limit laws for the distribution of X n , maximizing and minimizing asymptotically the 
value of the fourth moment E {X*) 4 , respectively, with other limit laws lying in between. 

A standard example where Gaussian limit law arises is the number of inversions in random 
permutations (or Kendall's r-statistic) 

p (^= n — i — • 

A straightforward calculation shows that the fourth normalized central moment has the form 

9(6n 2 + 15n + 16) 
~ 25n(n- l)(n + 1)' 

which implies the asymptotic normality by Theorem 1.1; see Feller (1945), Section 4 for more 
details and examples. 

On the other hand, a Bernoulli limit law results from the simple example 

1 4- z 2n 

It is then natural to ask to which limit laws other than normal and Bernoulli can the sequence 
of random variables X* converge. The simplest such example is the uniform distribution 

„ , , 1 + z + z 2 + • • • + z 2n 
y ! 2n + l 

or, more generally, the finite sum of uniform distributions 

1 + z H h z 2n i 



P ni +-+n k (z) — ][ 



2n j + 1 

Observe that the moment generating functions of the above three distributions have the 
following representations. 

- Normal distribution: e s2//2 ; 



- Bernoulli distribution (assuming ±1 with equal probability): 

e s + e~ s , A 4s 2 



cosh = n (i + {2k 4 _ s iy7r2 ) ^ 



- Uniform distribution (with zero mean and unit variance): 



JV3 V3lS ]£\ k 2 7T 2 J 
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Here we used the well-known expansions (see Titchmarsh (1975)) 

(4t5^ A sin s ( "\ 

1 - (2k - 1) v» J ' — = ni v i -^j- 

We show that these are indeed special cases of a more general representation theorem for the 
limit laws. 

Theorem 1.2. Le? {X n } Z?e a sequence of random variables whose probability generating 
functions are polynomials with only roots of modulus one. If the sequence {X*} converges 
to some limit distribution X, then the moment generating function of X is finite and has the 
infinite-product representation 

E(e Xs ) = e** 2 J] (l + f * 2 ) , (4) 
where q and the sequence {q k } are all non-negative numbers such that 

fcjsl 

The above examples show that q k = ^^l-i) 2 ^ or Bernoulli distribution and q k = for 
the uniform distribution. More examples will be discussed below. 

It remains open to characterize infinite -product representations of the form (4) that are 
themselves the moment generating functions of limit laws of root-unitary polynomials. On the 
other hand, many sufficient criteria for root-unitarity have been proposed in the literature; see, 
for example, the book Milovanovic et al. (1994) and the recent papers Lalin and Smyth (2012); 
Suzuki (2012) for more information and references. 

This paper is organized as follows. We first prove Theorem 1.1 in the next section when n 
is even, and then modify the proof to cover polynomials of odd degrees. Theorem 1.2 is then 
proved in Section 3. We then apply the results to many concrete examples from the literature: 
Section 4 for normal limit laws and Section 5 for non-normal laws. A class of polynomials 
with non-normal limit law is included in Appendix since the root-unitarity property has not yet 
been proved. 

2 Moments and the two extremal limit distributions 

For convenience, we begin by considering (general) polynomials of even degree with all their 
roots lying on the unit circle 

P2n(z) = P kZk ' 
0<fc<2n 

where p k ^ 0. To avoid triviality, we assume that not all p k s are zero. Observe that if \p\ = 1 
and P(p) = 0, then P(p) — 0. If p — 1, then its multiplicity must be even since all other roots 
can be grouped in pairs and are symmetric with respect to the real line. Thus our polynomials 
can be factored as 

where \pj\ = 1 for j = 1, . . . n. This factorization implies that root-unitary polynomials are 
always self-inversive. 
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Lemma 2.1. The coefficients of a root-unitary polynomial of even degree 2n are symmetric 
with respect to n, that is 

Pn-k=Pn+k (O^k^n). 

Proof. By replacing z by 1/z, we get 

P2n-kZ k = Z 2n P 2n (l/ Z )= Y[(l-Z Pj )(l-Zp]) 

0<fc<2n l<i<n 

= Yl (Z - Pj)(z -p]) = P 2n ( Z ) = ^ PkZ k - 
l<i<n 0<ifc<2n 

Taking the coefficients of z k on both sides, we obtain P2 n -k = Pk for ^ k ^ 2n, which 
proves the lemma. ■ 

2.1 Random variables, moments and cumulants 

Since the coefficients of P2n(z) are nonnegative, we can define a random variable X 2n by 

e(^) = ^44- 

For convenience, we write pj = since \pj\ = 1. Then 

— Pj)(z — pj ) = 1 — (pj + p~])z + Z 2 = \ — 2z COS 0j + z 2 . 

It follows that 

. y , -r-r 1 — 2,2 COS 0o + -2 2 

v ; 11 2(1- cos 0,-) 

Note that fa ^ for 1 < j < n since P 2n (l) > 0. 

It turns out that the mean values of such random variables are identically n. 

Lemma 2.2. For n > 1 

E(X 2n ) = n. (6) 

Proof. By (5), take derivative with respect to z and then substitute z = 1. ■ 

The relation (6) indeed holds more generally for self-inversive polynomials; see, for exam- 
ple, Sheil-Small (2002). 

Corollary 2.3. All odd central moments ofXm are zero 

E(X 2n -n) 2m+1 = (m = 0,l,...). 
Proof. This follows from the symmetry of the coefficients p^. ■ 
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For even moments, we look at the cumulants, which are defined as 
where K 2m +x{n) = 0. 

Lemma 2.4. The 2m-th cumulant K 2m (n) of X 2n is given by 

K2m{n) = (2m)! ^ — — — h m , k S n , k (m ^ 1), 

w/z<?r<? 2 2A: sinh 2/c (s/2) = ]>] m>fc h m , k s 2m , with h k)k = 1, and 

s n ,k ■= ?! 

1 — COS ( 

Proof. By (5), we have 

, 1 - 2e s cos + e 2s , / , (e s - l) 2 
log — — — = log e s + 



3 / 



k ' 



2(1- cos 0) °V 2(1- cos 4 

sinh 2 (s/2) 



= s + log 1 + 2- 

\ 1 — COS ( 

Thus 

,2/ 



logE(e^)^log(l + 2^^) 
£?< n V (1- cos 0,0 / 



which implies (7) by a direct expansion. 



(7) 



2.2 Variance and fourth central moment 

In particular, we obtain, from (7), 



„;:=V(* fc ) = *(„)= £ j-* (8) 



Lemma 2.5. 77ze variance satisfies the inequalities 



\ < < - 2 - (9) 



iVoo/ The lower bound follows from (8) and the inequality 1 — cos 0, ^ 2. The upper bound 
is also straightforward 



2r) \ - . 



which shows that the distance of the unit zeroes of P-miz) to the point 1 is always larger than 
c/n, where c > is an absolute constant. ■ 
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On the other hand, by the elementary inequalities 



2 t 2 

-t 2 — (t G [ — 7T, 7r]), (10) 



7T 2 ^2 

we have 

cr 2 7T 2 

2 < — ^ —■ 

We now turn to the fourth central moment. Define 



Lemma 2.6. (i) For n ^ 1, 



l^E — ^3- — <3. (12) 



W 

4 

— n \ 

Proof. By definition and by (7), 



X 2 „. n ^ 3 ^ Wn _^ a (J3) 



E = 3 + = 3 + - 3 ^n. (14) 



Now 

_ 2 1 2 + COS0j 1 1 x _ 2 



cos 6 a)' , 



proving the upper bound of (12). On the other hand, since l/a n ^ y/2/n (by (9)), we see that 
(13) also holds. It remains to prove the lower bound of (12), which results directly from the 
Cauchy-Schwarz inequality 



, Xo„ — n \ 2 / / Xo n — n x 4 x 
1 = E — I < E 



C n / \ \ cr 



By (10), we can replace the condition cu n — > by 



e^=°((e< 

!<.7<n \ l<.j'<n 



, 2 

-2 



On the other hand, K((X n — n) / a n ) 4 — > 3 is equivalent to n 4 {n) / n 2 (n) — >■ 0; the latter condition 
is in many cases easier to manipulate. 

Note that (12) proves (1) when n is even. 
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2.3 Estimates for the moment generating functions 

Lemma 2.7. Assume u n — > 0. For all s 6 C swc/z that \s\ ^ min{o„, w n 1//4 }/4, we /zave 

E(e (X2 "- n)s/CT ") = exp + ° (^- + w "l s l 4 )) ■ ( 15 ) 



Proof. By (5), 



Note that, by (8), 



2 1 

a„ ^ max 



l<i<n 1 — COS (f)j 

From the definition (1 1) of uo n , we also have 

1 1 

— max — - ^ u n , 

a* Kj<n (1 - COS0j) 2 



which means that 



Thus 



max — — ^ crly/uJ^. (16) 

l^jKn 1 - cos <pj 



e s/<r n _ 1+ 1) 



2(1 — cos cj)j) 



0"n 



Since cu n — >■ by assumption, the right-hand side is less than, say 2/3, for large enough n when 
\s\ remains bounded. Thus we can use the Taylor expansion of log(l + w) and obtain 

/ (pS/lTn _ 1 ^2 

log 1 + (e s/ ^ - 1) + 



2(1 — cos 4>j) 



s s 2 



+ 7^+0 —+ -r- — r + 



a n 20-^(1 — cos0j) V "^! — cos0j) o^(l — COS0-,) 2 o^(l — cos0j) 3 / 

By (16) 



|s| 6 A /cX 



0-6(1 - COS0J-) 3 ^ «7*(1 - COS0,-) 2 ' 

It follows, after summing over all j, that 

\ogE(e x ^) = t + - + O ( ^ + (\s\ 4 + \s\ 6 ^)u n ) . 

Now if \s\ ^ min{o n , w^ 1 ^ 4 }/4, then w^ 2 |s| 6 ^ w n |s| 4 /16, and this proves (15). ■ 
Lemma 2.8. For set, the inequality 

E ( e (Xan-nW<r») ^ gxp (| ^2*/^) (17) 

/lo/ds. 
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Proof. By (5) and the elementary inequality 1 + y ^ e y for real y, we obtain 



E<^)=ji(i + <,-i) +5 J^y 



Thus 



p n( Z -l)+^(«-l) 2 /2 



E(e (X2 "- n)s/CT ") *C exp ( n ^e s/<T " - 1 - ^-J + \o 2 n {e slcJn - if 

^ot 1 



and (17) follows from the inequality nja\ ^ 2. 



2.4 Normal limit law 

We now prove the second part of Theorem 1.1 in the case of polynomials of even degree, 
namely, {(X n — n)/o n } converges in distribution and with all moments to the standard normal 
distribution if and only if 

V — (7 J 

Proof. Consider first the sufficiency part. By (13), tu n — > 0, and we can apply the estimate 
(15), implying the convergence in distribution of (X 2n — n)/o n to jV (0, 1). 
On the other hand, by Lemma 2.8, 

^-2n — n \ S 3 2 p 2s/CTn 



o n ) (2m)! 

Taking s = 1, we conclude that all normalized central moments of X 2n are bounded above by 

„ f Xo n — 7l\ . , . 3„2/<r n 

E[— ^ (2m)!e2 e . 



Or. 



Thus we also have convergence of all moments. 

For the necessity, we see that if {(X 2n — n)/o n } converges in distribution to ^V(0, 1), then 
the the fact that the moments of (X 2n — n)/o n are all bounded implies that all the normal- 
ized central moments of X 2n converge to the moments of the standard normal distribution; in 
particular, the fourth normalized central moments converge to 3. ■ 
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2.5 Bernoulli limit law 

We now examine the case when the fourth moment converges to the smallest possible value, 
that is 



Note that 




If (18) holds, then by Chebyshev's inequality, we see that 

P { ^ 2n ~ n e ("I - e, -1 + e) [J(l - e, 1 + eU -> 1, 

for any e > 0. By symmetry of the random variable X 2n — n 

G (-1 - e, -1 + e)) = P (^-Z^_ g (1 — e, 1 + e)^ • 

We conclude that the distributions of (X 2n — n)j a n converge to a Bernoulli distribution that 
assumes the two values 1 and —1 with equal probability. 




2.6 Polynomials of odd degree 

To complete the proof of Theorem 1.1, we need to address the situation of odd-degree polyno- 
mials. 

Assume Q 2n -i(z) is a root-unitary polynomial of degree 2n — 1 with non-negative coeffi- 
cients. If we multiply it by the factor 1 + z, then the resulting polynomial 

P2n(z) = (1 + z)Q 2n -l{z) 

remains root-unitary with non-negative coefficients. This means that the moment generat- 
ing functions of the corresponding random variables E(e y2,I ~ lS ) := Q2n-i(e s ) / / <52n.-i(l) an d 
E(e X2nS ) := P 2 n(e s ) /-P2n(l) are connected by the identity 

E(e X2 " s ) = ^-E(e y2 - lS ). 

This leads to the relation 

X 2n = y 2n _! + B, (19) 
where B is independent of Y 2n -i and takes the values and 1 with equal probability. Thus 

E(Y 2n _ x ) = E(X 2n ) - ± = n - §, 

V(F 2 „_ 1 ) = V(X 2n ) -\ = al-\, (20) 
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and 

E (r 2n _x - E(F 2 „_ 1 )) 4 = E(X 2n - nf -lo* + ± (21) 

Thus we obtain 



2"n 16 



which, by (12), is bounded above by 



K-|) 2 V <? 2 J K-J) 2 A A(*A-i) 2 

Thus the fourth normalized central moment is bounded above by 3; the lower bound follows 
from the same Cauchy-Schwarz inequality used in the even-degree cases. 
On the other hand, since (again by (19)) 

X 2n - E(X 2n ) A vM^n-i) Yin-x ~ E(y 2n -i) B-\ 



V^X^J ^W^n) y/V(Y 2n -i) ^W^Y 
we have, by (20), 

X 2n - E{X 2n ) d F 2n _! - E(Y 2n _ 1 ) / „ / 1 



The last identity implies that both sides converge to the same limit law. 
Assume that the fourth central moment of Y 2n -i satisfies 

\ 4 

Y 2n -\ — E(l 2ra _ 



+ (22) 



Then, by (21), we obtain 
E 



X 2ra -E(X 2n ) \ / v(y 2 n-i) y f ^i-EM \ , _n 

Thus the left-hand side also tends to 3 and, consequently, X 2n is asymptotically normally dis- 
tributed. The asymptotic distribution of X 2n then implies, by (22), that of Y^n-i- 
The proof for the Bernoulli case is similar and is omitted. 



3 The infinite-product representation for general limit laws 

We first prove Theorem 1.2 in this section, and then mention some of its consequences. 
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3.1 Proof of Theorem 1.2 



It suffices to consider only the sequence of polynomials of even degree. The symmetry of distri- 
bution of the limit law X follows from the symmetry of coefficients of polynomials P 2n (z). The 
inequality (17) for the moment generating function of (X 2n — n)/a n implies that the moment 
generating function of the limit distribution X is also finite, and thus X is uniquely determined 
by its moments. This means that the sequence {(X 2n — n)/a n } converges in distribution to X 
as n — > oo if and only if 

E H^J^T E(X m ) (m > 0), 

as n — > oo. Thus the cumulant R m (n) of (X 2n — n)/a n of order m also converges to the 
cumulant of X of order m for m ^ 1. Note that R 2m +i( n ) = for m ^ and (see (7)) 

_ / \ _ 2m , N (2m)! ^ (-l) fc_1 . c 
K 2m {n) = o n K 2m {n) = ^ n mtk S nik . 



a 2m ^ 2 fe 



Since S„ k ^ o"^ fc , we the deduce that 



n,k 

m—l 



(2m)! m2 m 
for any fixed m. Now a n — > oo, we conclude that 



K>2m t K>2m(n) ( — l) m 1 S n m 

hm — — = hm ' . (24) 

/ C\ \ 1 O ryi 'Jim > ' 



(2m)! n->oo (2m)! m2 m n->oc cr 2m 
We now introduce the distribution function 



-n 1 <x 

t£(1— COS •) 



(T^l-COS^) 



with support in the unit interval. Then 



S n ,N I jv-l i;', i 



fl 



The fact that the left-hand side of the above expression has a limit (24) implies that the corre- 
sponding sequence of distribution functions F n (x) also converges weakly to some limit distri- 
bution function Fix). Therefore 



Tin 



lim ^ = / x^dF(x), 

n-»oo at 







which implies that the cumulants of the limit distribution X can be expressed as 

K 2m r R 2m {n) (-I)™' 1 S n , m (-l)™" 1 r 1 , 

r- = hm - — Vr = lim -^r— = / x dF(x). 

(2m)! n^oo (2m)! m2 m n->oo o^ 1 m2 m J v ' 
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It follows that 



E(e Xs ) = exp V ^s 2m 



/ ( _-\\m—l'lm /•! 



(2m)! 

m-l „2m /"l 



(25) 



/r-logqW/2) iF(x 

\Jo x 

Note that the distribution function F n (x) has no more than [l/e\ points of discontinuity in 
the interval [e, 1] if e > 0. Thus the weak limit F(x) of the sequence of F n (x) has the same 
property: F(x) has no more than [l/e\ points of discontinuity q k in the interval [e, 1], where 
qk is the limit of certain points of discontinuity of function F n (x). This means that F(x) is a 
distribution function of the form 



F{x) 



Q + E gfc <s Qk, if x ^ 0, 
0, ifa;<0, 



where q k > with J2k^i Ik = 1 — 9- Here g equals the jump of the function at zero. Thus 
f 1 log(l + xs 2 /2) , , g 2 v^i , ^ 

Substituting this expression into (25), we obtain (4). This completes the proof of Theorem 1.2. 



3.2 An alternative proof of Theorem 1.2 

A less elementary proof of Theorem 1 .2 relies on the Hadamard factorization theorem (see 
(Titchmarsh, 1975, Ch. 8); see also Newman (1974) for a similar context). Indeed, assume that 
(X 2n — n)/a n converges in distribution to some limit law X, then the inequality (17) implies 
that 

|E(e Xs )| ^ e 3 l s l 2 / 2 (seC). 

In other words, it is an entire function of order 2. Hadamard's factorization theorem then 
implies that such a function can be represented as an infinite product 

E(e Xs ) = e As2+Bs l[(l - -V/", 
p ^ 

where p ranges over all zeros of the function of the left-hand side. On the other hand, the fact 
that all zeroes of the functions E(e^ X2n_n ^ //CTn ) are symmetrically located on the imaginary line 
implies the same property for E(e Xs ). This yields 

E(e Xs ) = e As2+Bs H(l + S ^) : 
for some real sequence tk > 0. Now E(X) = implies that B = 0. Also E(X 2 ) = 1 leads to 

Denoting by q = 2 A and qk = 2/t 2 ,, we obtain the representation (4). 
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3.3 Implications of the infinite-product factorization 

By (4), 

C_1 \m-l 



m2 r 

3>l 



This yields the sign-alternating property for the sequence {/t2m}- 

Corollary 3.1. If X is not the normal law, then all even cumulants are non-zero and have 
alternating signs 

(-l) m -^ 2m >0 (m^l). 



Corollary 3.2. 

Proof. By (4), 



1 *C E(X 4 ) ^ 3. 



E(X 4 ) = 3 l-]Tg 2 , (26) 
which implies the upper bound; the lower bound follows directly from Cauchy-Schwarz in- 



equality 1 = E(X 2 ) ^ v/E(X 4 ). ■ 

Corollary 3.3. The standard normal distribution is the only distribution for which the fourth 
moment reaches the maximum value 3 in the class of distributions that are the limits of random 
variables whose probability generating functions are root-unitary polynomials; similarly, the 
Bernoulli distribution assuming ± 1 with probability 1/2 each is the only distribution whose 
fourth moment reaches the minimum value 1 in the same class of distributions. 

Proof. Note that the standard normal law corresponds to the choices q = 1 and qj = 0, the first 
part of the corollary follows then from (26). 

For the lower bound, assume that Y is a symmetric distribution such that E(F) = and 
E(F 2 ) = E(F 4 ) = 1. Then 

Y(Y 2 ) = E(Y 2 - l) 2 = E(F 4 - 2F 2 + 1) = 0. 

This means that Y can only assume two values F(Y E { — 1, 1}) = 1. The symmetry of Y now 
implies that Y assumes the values 1 and —1 with equal probabilities. ■ 

Remark 3.4. The uniqueness of the standard normal and Bernoulli laws also implies that a 
sequence of random variables {X n } converges to normal or Bernoulli if and only if its fourth 
normalized central moment converges to 3 or to 1, respectively. This provides an alternative 
proof of the last two statements of Theorem 1.1. 

4 Applications. I. Normal limit law 

We consider in this section applications of our results in the situations when the limit law is 
normal. 
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4.1 A simple framework 

Our starting point is the polynomials of the form 

P (z) - (l-^Xl -**)••• (!-**") (21) 
A } (1 - z^){l - #*) ■ ■ ■ (1 - z*») ' K ] 

where aj, bj are non-negative integers that may depend themselves on N and 



£ ft 



n := > 



We assume that P n (z) has only nonnegative coefficients. Such a simple form arises in a large 
number of diverse contexts, some of which will be examined below. In particular, it was studied 
in the recent paper Chen et al. (2008). 

We now consider a sequence of random variables X n defined by 



E(z 



X n \ _ Pn{z) 



Pn{l) 



We have 



where 



pji) " exp U- m \ s 



i m 

— B ™ ( b T~ a T) ( m ^ X ) 



m 



the B m 's being the Bernoulli numbers. Note that £>2 m +i = for m ^ 1. 
An application of Theorem 1.1 yields the following result. 

Theorem 4.1. 77ze sequence of the random variables (X n — K(X n ))/ a/V(X„) converges to 
the standard normal distribution if and only if the following cumulant condition holds 

v K N,4 144 Ei^ K ^^ 4_a i) n 

hm = — — hm — — k = 0. (28) 



\Zul^N{°j a j) 



The cumulant condition largely simplifies the sufficient condition given by Chen et al. 
(2008), where they require the convergence of all cumulants 

— - > (m > 2), 

following the proof used by Sachkov (1997). See also Janson (1988) for a related framework. 



4.2 Applications of Theorem 4.1 

Theorem 4.1 can be applied to a large number of examples. Many other examples related to 
Poincare polynomials, rank statistics, and integer partitions can be found in the literature; see, 
for example, Akyildiz (2004); Andrews (1976); van de Wiel et al. (1999) and the references 
therein. 
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Inversions in permutations The generating polynomial for the number of inversions in a 
permutation of n elements (or Kendall's r statistic) is given by 

TT — ■ 

H l-z 

In this case, the cumulant condition (28) has the form 

(El^nO' 2 - 1)) 

Thus the number of inversions in random permutations is asymptotically normally distributed; 
see Feller (1945), Sachkov (1997); see also Cronholm and Revusky (1965); Louchard and 
Prodinger (2003); Margolius (2001). 

Number of inversions in Stirling permutations In this case, we have the polynomial (see 
Park (1994)) 

„ 1 _ 7 r+(j~l)r 2 

n rzp (rsi) ' 



and the cumulant condition (28) is of order 

■ + ?r 2 ) 4 - 1) 

0{n~ x ). 



K 



n,2 



(Zo^<n((r + jr*T - 1) 



Consequently, the number of inversions in random Stirling permutations is asymptotically nor- 
mally distributed. 

Gaussian polynomials The generating function for the number p(n,m,j) of partitions of 
integer j into at most m parts, each ^ n, is given by (see e.g. Andrews (1976)) 

Then the cumulant condition has the form 



(£l<i<n(( m + ~ 3 



2 

2\ ' 



m n 



This means that the coefficients of Gaussian polynomials are normally distributed if both 
n,m — > oo; see Mann and Whitney (1947); Takacs (1986). More examples can be found 
in Andrews (1976). 
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Mahonian statistics In this case the polynomials are equal to the general g-multinomial co- 
efficients (see Canfield et al. (201 1) and Canfield et al. (2012)) 



Pjz) 



where n = £ 2 <fc<m a & Ei<j<fc a r By symmetry, we can assume that a\ > - • • ^ a m . Then 
the cumulant condition (28) becomes 

Sl<i<ai+-+a m * _ El<j<m El<i<aj * /4( a l H + Q-m) ~ Elsysgm fi{ a j) 



(Y,<-< . . I 2 - El*-* ^ (/2(fll + ---+a m )-El^m/2W) 

\ ^— 'l^J^OiH him ZL-zlsgjigm /lsgisga., / 



2 • 



where / 2 (x) = (2x 3 + 3x 2 + x)/6 and /^(x) = (6x 5 + 15x 4 + 10x 3 — x)/30. By induction, 
(ai + • • • + a m ) k — a\ — ■ ■ ■ — a k m is nonnegative and is nondecreasing ink ^ 1. Thus the 
right-hand side is bounded above by 



5 



9-31 (di H + a m ) 5 — a\ a: 

30 ((ai H h a m ) 3 - a\ a%) 



2 



() , AH hflm \ _ q( a l H 1" a r 



J2ltZi<j^m a i a 3 ) \ a l( a 2 + Ct3H Ha 

oi—^ + 1 

a 2 + a 3 H h a m ai 



where we use the estimates 



(ai H h a m f -a\ a m x ( a i H 1- X] 

(ai + • • • + a m f -a\ a 5 m x (a x + • • • + a m ) 3 



Thus we arrive at the same conditions as those in Canfield et al. (201 1) 

a\ — > oo and a 2 + a 3 + ■ ■ ■ + a m — > oo, 
for the asymptotic normality of the coefficients of P n (z) when a\ ^ a 2 ^ • • • ^ a m . 

Generalized g-Catalan numbers The generating function has the form 

and the cumulant condition (28) also holds 

£ W „(((ro - 1)" + ^ - J 4 ) ^ E 2<3 -<n( 2mw ) 4 = q 



(E w „(((m - l)n + j) 2 - J 2 )) (E Wn (m - 1) 2 ^ 2 



which means that the generalized g-Catalan numbers are asymptotically normally distributed, 
uniformly for all m ^ 2. This result was previously proved by Chen et al. (2008). 
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Sums of uniform discrete distributions Let X n be the sum of iV independent, integer- valued 
random variables 

X n := J\ + J2 + ■ ■ ■ + Jni 
where J k is a uniform distribution on the set {0, 1, 2, . . . , d k — 1} with d k ^ 2, and n = 
Yli<j<N^3 ~ !)■ Then the corresponding probability generating function E,(z Xn ) is equal, up 
to a normalizing constant, to 

nV 



Pn{z) = 



z dj 



which means that X n is asymptotically normal if and only if 

•2 



Since by our assumption dj ^ 2, we have gL, — 1 x gL,- and thus we can simplify our necessary 
and sufficient condition for asymptotic normality as 

(4 + 4 + ... +d ir- ■* < w -> °°>- <29) 

Note that dj here can depend on N. The continuous version of this problem with J k being 
uniformly distributed on the intervals [0, dj] was considered in Olds (1952). The corresponding 
necessary and sufficient condition obtained in this paper was 

maxi^jvcL 
. -> 

y/dl + d 2 2 + --- + d% 

which is equivalent to condition (29). 



Number of inversions in bimodal permutations A permutation o = (sx, s 2 , . . . , s n ) of 

n numbers 1,2,3, ... ,n is said to be of a shape (i, k — j,j, Z) if the first i numbers in the 
permutation are decreasing si > s 2 > ■ ■ ■ > s», the next k — j numbers are increasing Sj + i > 
s 2 > • • • > s i+k _j, then followed by j increasing and Z decreasing numbers. Assume that a is 
chosen with equal probability among all permutations of shape (i, k — j, j, I). Then its number 
of inversions becomes a random variable I n = I n (i,k — j,j, I). The probability generating 
function of I n is, up to some constant, of the form (see Bohm and Katzenbeisser (2005)) 

The random variables I n are asymptotically normally distributed if 

+ ») 4 - » A ) + ELi((fc + j + ^) 4 ~ ^ 4 ) + ~ i + ^) 4 ~ ^ 4 ) , 0; 

+ - ^ + EU((* + * + ") 2 - "*) + £l=i((* - J + ") 2 - ^ 2 )) 2 
which is equivalent to 

ik(k + i) 3 + l(k + i)(k + i + I) 3 + j(k 4 - j 4 ) ^ 
(iife(A; + i) + Z(Jfe + z)(Jfe + i + Z) + j(^ 2 - j 2 )) 2 
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If we assume that the parameters k, I are proportionate to some parameter t, that is i = 
[at\ , j = [/3t\ , k = [_7tJ , I = [St\ , where a, /3, 7, 5 > and a + 7 + 5 = 1, then the above 
condition is satisfied and as a consequence I n is asymptotically normally distributed as t — > 00. 
This fact has been proved in Bohm and Katzenbeisser (2005) by the method of moments. 

Rank statistics Many test statistics based on ranks lead to explicit generating functions that 
are of the form (27), and thus the corresponding limit distribution can be dealt with by the tools 
we established. In particular, we have the following correspondence between test statistics and 
combinatorial structures; see van de Wiel et al. (1999) for more information. 



Kendall's r 


Inversions in permutations 


Mann- Whitney test 


Gaussian polynomials 


Jonckheere-Terpstra test 


Mahonian statistics 



On the other hand, the Wilcoxon signed rank test (see Wilcoxon (1947)) leads to the prob- 
ability generating function of the form 

which admits a straightforward generalization to (see van de Wiel et al. (1999) for details) 

n l + Z*i_ 
2 ' 

where the a/s can be any real numbers. When they are all nonnegative integers, we see, by 
(28), that the associated random variables are asymptotically normally distributed if and only 
if 

4 + • • • + a n n 

K + --- + 2 ' 

as n — > 00. In particular, this applies to Wilcoxon's test (a, = j) and to Policello and 
Hettmansperger's test (aj = min{2j, n + 1}; Policello and Hettmansperger (1976)). 

4.3 T\iran-Fejer polynomials 

The class of polynomials we consider here (see (30) below) is of interest for several reasons. 
First, they lead to asymptotically normally distributed random variables but do not have the 
finite-product form (27). Second, they provide natural examples with non-normal limit laws 
when the second parameter varies. Finally, they have a concrete interpretation in terms of the 
partitioning cost of some variants of quicksort. 

Fejer (1937) studied the Cesaro summation of the geometric series defined by 

F n , k {z) := F o,k-i{z) {k > 1), 

with 

F n , (z) := z ^ 



19 



and Turan (1949) proved that all F^l(z) are root-unitary for ^ k ^ n. We characterize all 

possible limit laws for the random variables defined via the coefficients of F^ k \z) for ^ k ^ 
n. 

By the relation 

F n>k (z) = [W n \- — — r, 

(1 — w) k+L (l — zw) 

where [w n ]f(w) denotes the coefficient of w n in the Taylor expansion of f(w), we have 



F w U) = \w n ] 1 ^ 

*n,kK Z ) l W \l- W )k+l (l- zw ) 



fe+1 



* e crxv) 



Normalizing this polynomial, we obtain 

P nM--= £ ^ (30) 

Q<j<n-fc V 2fc+l / 

which gives rise to a sequence of probability generating functions of random variables, say 
Z nk . Note that 



3\ n-l-j\ 

Z fn-k-x^Z) - 7 n \ 

fc<i<n-fc-l V2fe+1/ 



which arises in the analysis of quicksort using the median of 2k + 1 elements; see Sedgewick 
(1980); Chern et al. (2002) or Appendix. 

Lemma 4.2. For m ^ 

0<^<m V 2fc+l / 

where S(m, £) denotes the Stirling numbers of the second kind. In particular, 

E(Z n , fc ) = — anJ V(Z n , fc ) = ^ 4( A fc + 3) ^. (32) 

Proof. By (30), the relation 

and the combinatorial identity 



j + k\ fn — j\ f j\ fk+£\fn + k + l 

k )\2k + £ + l 



(easily proved by convolution), we deduce (31). 
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Theorem 4.3. The random variables Z n ^ k are asymptotically normally distributed if and only 
if both k and n — k tend to infinity. IfO ^ k = 0(1), then the limit law is a Beta distribution 

-^-^Beta(Jfe,ife). (33) 
n 

Ifl ^ I := n — k = 0(1), then the limit law is a binomial distribution 

Z ntk -A Binom(£; \). 

Proof. By (32), the variance tends to infinity if and only if n — k — > oo (0 ^ k ^ n). Also we 
obtain, by (31), 

E (Z n;k - ^) 4 _ _ 2(3n 2 + Qn + k 2 + Ak + 6) = ( n 



Y(Z n , k ) 2 (n - k)(n + k + 2)(2k + 5) \k(n - k) 

The asymptotic normality then follows. We can indeed obtain a local limit theorem by straight- 
forward calculations from (30). 

When k = 0(1), we have, by (30), 

E TO ^ (fc + m)!(2fc + l)! 

n m k\(2k + m + l)\ ^ n ^ v )> 

implying that the moment generating function of the limit law satisfies 

(2Jfe + l)!y> (k + m)\ n , m _ (2k + l)\ f 1 

m>0 



ZhS . (2k + l)\^ (k + m)\ m (2k + l)\ f 1 kn . kxS j 

E(e ZkS ) = - — -—- > — — ^ s m = v —- 1 I x k (l - x) k e xs dx, 

K ' k\ ^ m\(2k + m+ 1)\ k\k\ J 



a Beta distribution. Note that we can express the moment generating function in terms of Bessel 
functions as 

E(e^-W) = (^y k " /2 T(k + l)J k+h (is/2) 

m \ ^ k+1 / 2 ,jJ 

where J a denotes the Bessel function and the Ca,/ S denote the positive zeros of J a (z) arranged 
in increasing order. By considering 2(Z k — l/2)\/2k + 3, we obtain (4) with qj = 2(2k + 

3)/ d+i/2,j- 

On the other hand, when £ := n — k = 0(1), we have, by (31), 



P n .k(z) ->■ 



1 + z 



2 

a binomial distribution. Note that we have the factorization 

As 2 

(2 J - 1) 2 7T 2 £ 
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5 Applications II. Non-normal limit laws 



In addition to the extremal cases of the Turan-Fejer polynomials, we consider in this section 
more root-unitary polynomials whose coefficients have a limit distribution that is not Gaussian. 
A class of polynomials exhibiting a similar non-Gaussian behavior is included in Appendix 
because the proof that they are root-unitary is still missing. 



5.1 Reimer's polynomials 

In the course of investigating the remainder theory of finite difference, Reimer (1969) proved, 
as a side result, that the polynomials 



Rn, m (y) ■= E ("V / 



i+i 

\t(t- l)---(t-n- l)\ m dt. 



have only unit roots. We consider the distribution of the coefficients of R n>m (y). 

For simplicity, we consider only m = 1 and write R„ = R rh \. Define the random variables 
X n by 

E (^») := ^M- 

1 ' Rn(l) 

Let 

log(l - z) 

These numbers are (up to sign) known under the name of Cauchy numbers; see (Comtet, 1974, 
pp. 293-294). See also the recent paper Kowalenko (2010) for a detailed study of these num- 
bers. 



Lemma 5.1. For n ^ 1 



E(y*») = 12 E ( n ) V n ~ j i} ~ v)'>W ( 35 ) 



Proof. We have 

ri+i 



Rn(y) = E ( n ) (-V n+1+j y j r t(t-i)--.(t-n-i)dt 

^K-irEQi-wfCy* 



n 2)!(-l)" +L |.:" + 1 > I )(-!)'//■' /' +1 (l + ^dt 



+ 2)!( _ l) n + l [ ,n +2] £ ^V-1)V / 

+ 2)!(-l)"+ 1 [* B+1 ] E ("VWr^ x 



0<.?<n 

(n^)!^] ' 1 -^" . (36) 
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In particular 

R n (l) = (n + 2)\[z] ' -<" + 2 >' 



'log(l-z) 12 

and (35) follows. ■ 
Note that A = — 1 and 

_A< 
k- 

0^j<k 

All A k s are positive except A . 
Lemma 5.2. The moments of X n satisfy 

E(X™)= ^ A k S(m, k)n{n - 1) • • • (n - Jfe + 1) (to > 0), (37) 



A k = - y — ^ — (A; ^ i). 

A^. k + l-j ^ 



(38) 



O^i^k 

In particular, 

Ti Ti 

E(X n ) = -, V{X n ) = —(4n + U). 

Proof. By taking m-th derivative with respect to y and then substituting y = 1 in (36), we 
obtain 

,9™ n — d — 7)nY 

E(X n (X n -l)---(X n -m+l)) = 12[^ +:l1 



'<9?/ m log(l-*) 



y=l 



= I2[^ m+1 ]-^— J— n (n - 1) • • • (n - m + 1), 
log(l - z) 

which yields (37) since 

E(X™) = S(m,k)E(X n (X n -l)---(X n -k + l)). 



Theorem 5.3. The sequence of random variables {X n / n} converges in distribution to X whose 
m-th moment equals A m (defined in (38) J. 

Proof. By (37), E(X™) ~ A m n m . Since A k = 0(l/k), we see that A m = 0(2 m ), implying 
that such a moment sequence determines uniquely a distribution. ■ 
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The limit law has the moment generating function 



27TZS Viogl 1 + f ) S 2 

where the integration is taken along some Hankel contour; see (Flajolet and Sedgewick, 
2009, p. 745). 

When m ^ 2, we can apply the same arguments but the technicalities are more involved. 

5.2 Chung-Feller's arcsine law 

The classical Chung-Feller theorem states that the number of positive terms W n of the sums 
S n = X\ + • • • + X n , where Xj takes ±1 with probability 1/2 each, has the probability 

«w =*>=(?) ft (*=o,....»). 

The limit distribution is an arcsine law (see (Feller, 1968, §111.4)) 

— ■ - — W, where F(W < x) = — arcsin ^fx. 
n n 

The corresponding probability generating function is a polynomial with only unit roots. 
Indeed, following the same proof as in Turan (1949), we can show that K(z Wn ) is connected to 
Legendre polynomials by the relation 

E(z Wn ) = \v m 1 



y/(l-v)(l-zv) 

z l/2 + z -l/2 



= z n ' 2 Legendre r 

so that the root-unitarity of the left-hand side follows from the property that Legendre polyno- 
mials have only real roots over the interval [—1,1]. Note that the moment generating function 
of the arcsin law with zero mean and unit variance is given by the Bessel function 

2k\ (s/vzy 



E(e (W-l/2)./V5) = e -V~2s + (T) 



= j^is) = n ( i + % ), 

where the Co,/ s are the positive zeros of Jq(z). So we have (4) with q = and qj = 
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In a more general manner, from the Gegenbauer polynomials, one can also define the ran- 
dom variables W n by 

E(z w ») = ,„_,,!, ?\ b"i 1 



(2«+«-l) 1 J (1 - v)~ a (l - Zv)~ a 
/a+j-l\ fa+n-j-l\ 

= E J 'A-V * j («>o), 

for which all coefficients are positive and E(^ w/n ) has only unit roots. The limit law W a can be 
derived as in the bounded case of the Turan-Fejer polynomials 

E^-V*).) = f™\ ° +I/2 r(a+l/2)J„_ 1/2 (i.s/2) 

=n(i+ 4C , 

Note that the random variable 2y/2a + l(W a — 1/2) has variance one. 
For other potential examples, see (Johnson et al., 2005, Chapter 6). 

5.3 Uniform distribution 



s 2 



The literature abounds with criteria for the root-unitarity of polynomials. Among these, Lakatos 



(2002) proved that a complex polynomial P(z) := J2 Q<k<n a k z k with a k = a„_ fc is root-unitary 



see also Schinzel (2005). In particular, if the coefficients of P(z) are close to a constant, then all 
its roots lie on the unit circle. For example, let Ej = j![V]l/ cosh(^) denote Euler's numbers; 
then the polynomial 



PJz) 



cos(a/u7) cos(v / mJz) 



is root-unitary (see Lalin and Rogers (201 1)) with non-negative coefficients. See also Lalin and 
Smyth (2012) for more information and other root-unitary polynomials. Observe that 

as j, n — j — )• oo. Thus we can show that the random variables associated with the coefficients 
of P„(z) will be close to uniform, and the limit law is also uniform. Details are omitted here. 
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Appendix. A class of mixtures of hypergeometric distributions 

Yet another class of polynomials with a similar nature to those of Turan-Fejer arises from the 
analysis of the partitioning stages of quicksort and defined as follows (see Chern et al. (2002)). 
Consider the random variable Y n defined by 



where r ^ 1 and J2o<j<rPj = 1 is a given known distribution. Many concrete examples 
are discussed in Chern et al. (2002). Let P(z) : = J2o<j<rPj z:> • Assume throughout that 
Pj = Pr-i-j for 1 ^ j < r. Numerical evidence suggested that the probability generating 
function 



10:97-122. 



F(Y n = k) = Pi 



0^j<r 
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are root-unitary for many natural choices of {pj}, but it is unclear for which class of poly- 
nomials P(z) will the polynomials be root-unitary 2 . The results below do not depend on the 
root-unitarity of E,(z Yn ). 

Assuming from now on pj = p r -i-j for ^ j < r, we examine the moment structure of 
Y n . Note that this assumption implies that J2o<j< r JPj = ( r ~~ l)/2- 
Lemma 5.4. The m-th moment ofY n is given explicitly by 

(n+h\ 

m?) = E (39) 

Here 

v m ,h:=(-l) m+h S(m+l,h + l) s(/i + M + 1K 

where the s(m, h) denote the signless Stirling numbers of the first kind, and ii£ := ^2 ^j <r Pjj e - 
In particular, E(Y n ) = (n — l)/2 and /7ze variance satisfies 

_ (4tt 2 - r 2 + 3r)n 2 + 2(6vr 2 - 2r 2 + 3r - l)n + 8vr 2 - 3r 2 + 3r - 2 
( n) ~ 4 (r + l)(r + 2) ' 

Note that 7r 2 — 4r 2 + 3r is always positive because 

r(r ~ 3) , (r-1) 2 

~ T~ < = — ^ < vr 2 . 

Thus the limit law of Y r , is never Gaussian for finite r. 



Proof. By definition, 



fk\ /n—l—k\ 

ecc) = E Pi E fcm - ff" 3 • ( 4 °) 



0^j<r 

The, using the relation 

k m = J2 (-l) m+h S(m + l,h + l) (k + !)••• (k + m), 



we obtain 



n — 1 — A; 
j7 V r- 1 - j 



E (-ir + ^(- + i,h + 1) j> + 1) • • • (* + h) ( k ) ( n \ _ k \ 

^ (-i)^s(m +!,/»+ do- + 1) . . . o + ft) £ (* + *) (; : 1 : j) 

E (-l) m+h S(m + 1, + l)(j + 1) ■ • ■ (j + ri) + f) . 



2 The problem can be formulated by asking for which class of polynomials P(z) = J2o<j<rPj z '' w ^ tne 
polynomials 



\w n ] 1 P I 1 

L J (l- w){l- zw) \{l-w)(l-zw) 



have only unit roots? 
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Now, by substituting the expression 

(j + l)...( J + h)= s(h + l,£+l)j\ 

We then obtain (39). ■ 

Theorem 5.5. The sequence of random variables {Y n /n} converges in distribution to a limit 
law Y whose moment generating function satisfies 

E{e Ys ) = r Pj^ ~ X ) / e xs x r - x - j {l - x) j dx. (41) 



0^j<r 

Proof Indeed, by (39), 



J 



n m ( r + m)\ ^ 



so that 

where the moment generating function of X is given by 

-,m 



_ s m ^ r(r — j + m)r(r + 1) 
" m \ Z_> p i Y(r + m + l)T(r - j) 

= r E S £ Pi C 7 _ a; ) J ' a;r " 1 " J+,nda; ' 

m^O ' (Kj<r ^ / •'O 

which proves (41). The justification of the unique characterization of this limit law is straight- 
forward. ■ 

Examples. 

• The median of (2k + 1) elements: This corresponds to the case when r = 2k + 1 and 



Pk — 1. We then obtain 



E{e Ys ) = (2 ^, 1)! I' e xs x k {l-x) k dx, 



k\k\ 

a Beta distribution (with a Bessel-type infinite-product representation); see (33) and (34). 

• Uniform distributon: In this case, Pj = 1/r, 1 ^ j < r. This is algorithmically uninter- 
esting, but has the limit moment generating function (e s — 

• The ninther (the median of three medians, each being the median of three elements): This 
is the case when 

= {o,o,o, 0,0,0}. 

We have a mixture of Beta distributions for the limit law. 
Many sophisticated cases can be found in Chern et al. (2002). 
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